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ABSTRACT 



It is usually assumed, in classical statistical mechanics, that the tem- 
perature should coincide, apart from a suitable constant factor, with the 
mean kinetic energy of the particles. We show that this is not the case for 
Fermi-Pasta-Ulam systems, in conditions in which energy equipartition 
between the modes is not attained. We find that the temperature should 
be rather identified with the mean value of the energy of the low frequency 
modes. 
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1 Introduction. 

The Fermi-Pasta-Ulam system consists of a chain of N nonlinear oscillators 
with certain given boundary conditions, tipically fixed ends. It is well 
known (see jT| and |2|) that, for energies below a certain threshold Ec, if 
the energy is initially given to a few low frequency modes, equipartition 
of energy among the modes is eventually attained only after an extremely 
long time, while at intermediate times a kind of metaequilibrium state is 
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attained, in which the energy is shared essentially within a packet of low 
frequency modes. 

An interesting and much discussed problem, is whether the specific 
energy threshold Ec/N vanishes or not in the limit of an infinite number 
of particles. Here we leave this problem aside: we will suppose for example 
that the number N of particles be fixed, so that the threshold certainly 
exists. We address instead the problem raised by the fact that below the 
threshold one meets, as in the theory of glasses, with time scales to thermal 
equilibrium which are very long, even longer than any fixed observational 
time scale. Does this lack of thermalization have any consequences on the 
relevant thermodynamical quantities? or even, is it possible to correctly 
(i.e. uniquely) define the quantities of interest? In other words, is it still 
possible to have a thermodynamics below the threshold? 

In the literature, the discussion is usually focused on the specific heat, 
because from heuristic arguments it is suggested that to less chaotic mo- 
tions there correspond smaller specific heats, with eventually zero specific 
heat for totally ordered motions (i.e. for integrable systems). Thus it is ex- 
pected that by lowering the energy below threshold the specific heat should 
diminish; such a property, in turn, should be considered as a good indicator 
of the weakening of chaos. 

The papers (2j and 14), which aim at evaluating the specific heat of 
Fermi-Pasta-Ulam systems below threshold by numerical computations, 
reach two opposite conclusions: in the value of the specific heat remains 
constant (as would follow from the equipartition principle) even below the 
threshold, while in 0] the specific heat indeed begins to fall down, be- 
low the threshold, and finally vanishes as temperature goes to zero. This 
striking difference is apparently due to the different methods used in the 
two papers in estimating the specific heat. Actually, in both papers the 
Fermi-Pasta-Ulam system is kept isolated (fixed ends), so that a direct 
measurement of the specific heat is precluded (because a direct measure- 
ment requires at least one heat bath). The specific heat is then estimated 
from the fluctuations of energy of a subsystem through the well known 
relation between specific heat and mean square deviation of energy, which 
holds in the canonical ensemble. The two papers differ in the choice of the 
subsytem: in the paper "7 , one considers the energy fluctuations of a small 
piece of the total chain, while in W one considers the energy fluctuations 
of a small packet of nearby modes. As the energy of each mode remains 
nearly constant below the threshold while the energy of a piece of chains 
still present large fluctuations, this indeed explains why the two papers 
reach opposite conclusions. Now, at most only one of the conclusions can 
be correct, if a right conclusion does exist at all; indeed it is not clear 
whether the specific heat can be defined in an unambiguous way below the 
threshold. 
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A different approach was proposed in the paper 5 . In short, the basic 
remark is the foUowing one. As the above mentioned relation between spe- 
cific heat and energy fluctuations is obtained from the equilibrium Gibbs 
ensemble, then its validity below the threshold is in doubt just because, 
up to the considered times, the system has not yet thermalized. So, it is 
argued that in measuring the specific heat one should revert to the conven- 
tional method which makes use of a heat bath at a given temperature T 
coupled to the Fermi-Pasta-Ulam system, with the corresponding familiar 
calorimetric expression for the specific heat. Namely, the energy exchange 
AQ between the bath and the Fermi-Pasta-Ulam system is measured when 
the temperature is varied by AT, then the ratio AQ/AT is computed, and 
(in principle) the limit is taken for vanishing AT. However, even with 
such a method, one is still confronted with a delicate problem, because the 
amount of exchanged energy AQ does depend on how much time one has 
waited in making the measurement (this is the so-called waiting time of 
the theory of glasses). The curve predicted by equipartition is recovered 
for infinitely long waiting times, while for finite waiting times the specific 
heat is expected to vanish at sufficiently low temperatures. 

The question is thus: does there exist a natural way to choose a definite 
waiting time? Another problem then arises, due to the fact that, in the 
ratio AQ/AT defining the specific heat, one should insert in the denom- 
inator the variation of temperature of the Fermi-Pasta-Ulam system and 
not that of the bath. The question is then whether the temperature of 
the Fermi-Pasta-Ulam system is the same as that of the heat bath. The 
very fact that the quantity AQ depends on the waiting time actually shows 
that this is not the case, just because the equality of the two temperatures 
would imply AQ = 0. On the other hand, if one were able to identify 
the temperature of the Fermi-Pasta-Ulam system, then the question of 
the waiting time would have a quick answer: one should wait until the 
heat bath temperature and that of the Fermi-Pasta-Ulam system have be- 
come equal, and only at that time should one measure the corresponding 
energy exchange. So the possibility of having available well defined ther- 
modynamic quantities on short time scales is based on the possibility of 
providing a good notion of temperature for the FPU system before com- 
plete equipartition be achieved. From this point of view, the identification 
usually made of the temperature of the Fermi-Pasta-Ulam system with the 
mean kinetic energy of its particles is not the correct one, because in such a 
case the temperatures of the two systems (Fermi-Pasta-Ulam system and 
heat bath) would remain different for extremely long times. 

The identification of the mean kinetic energy with temperature is so 
deeply rooted in our minds, that the existence of another quantity playing 
that role seems hardly conceivable. The aim of this paper is to show instead 
that this is possible. 
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In Section 2 we give a preliminary discussion of the zcroth law for states 
of metaequilibrium such as those of the Fermi-Pasta-Ulam system below 
threshold, in Section 3 we describe the model we employ for measuring of 
the temperature of the Fermi-Pasta-Ulam system through heat baths by 
numerical computations, and the numerical results are given in Section 4. 

2 Zeroth law and temperature in states of 
metaequilibrium. 

One of the basic features of thermal equilibrium is the so-called zeroth law, 
which essentially amounts to the transitivity of the equilibrium. From this 
follows (see 0) that for any system there exists a function of its macro- 
scopic (the so called empirical temperature) which has the same value for 
bodies in equilibrium. 

However, it is not granted that, for a given macroscopic state, the equi- 
librium is unique if some of the internal degrees of freedom are dynamically 
frozen. We are thinking typically of the case of polyatomic molecules (see 
[S]) for which it is known that the exchanges of energy between the center 
of mass and the internal degrees of freedom are so slow that the number 
of effective degrees of freedom depends on the time of observation. This is 
actually the general situation that occurs in states of metaequilibrium. 

We thus address the problem whether it is possible to have a zeroth 
principle, and so also an empirical temperature, in situations of metaequi- 
librium, in which the physical quantities are changing only on a very long 
time scale. So, if we put our Fermi-Pasta-Ulam system in heat contact 
with another body, and observe that at first there is a rapid relaxation to 
a certain state, while a later evolution to a final equilibrium would take 
place on a time scale much longer than our observational scale, we can 
think of our system as if it were equilibrated. Obviously, one is not granted 
that in such a situation the zeroth law, i.e. the transitivity of this metae- 
quilibrium state, holds. But, if this is the case, an empirical temperature 
can be defined. In other terms, if we put the Fermi-Pasta-Ulam system in 
contact with a thermometric body which, after a short transient, appears 
to have reached a temperature T (not evolving on our time scale), and if 
later, after having put the Fermi-Pasta-Ulam system in contact with an- 
other body at the same temperature T, nothing seems to occur (i.e. there 
is no exchange of energy in mean between the bodies), then the metastable 
state reached does have the transitive property, and we are authorized to 
assign to the Fermi-Pasta-Ulam system the temperature T reached by the 
thermometric body. 

An equivalent arrangement, which we have actually implemented in 
our numerical simulations to be described below, is the following one: the 
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Fermi-Pasta-Ulam system is put at the same time in contact both with a 
heat bath and with a thermometric body (obviously, with no direct connec- 
tion between the two external bodies). In such a situation, the transitive 
property reduces to the property that in a short time the thermometer 
attains the same temperature of the heath bath; the subsequent evolu- 
tion to a final equilibrium should take place later, at a much smaller rate. 
Then, by definition, the temperature of the Fermi-Pasta-Ulam system in 
the metaequilibrium state is the one reached by the thermometer after the 
short-time relaxation. 

This definition can appear satisfactory from an operative point of view. 
However, as it stands, it still lacks a clear connection with the properties of 
the Fermi-Pasta-Ulam system itself. Indeed there remains the problem of 
understanding, how the zeroth law can hold even if the Fermi-Pasta-Ulam 
system did not yet thermalize. In this connection, we make reference to a 
known phenomenon [2| concerning the isolated Fermi-Pasta-Ulam system, 
namely the fact that, below threshold, for initial excitations of the low fre- 
quency modes the energy turns out to remain confined to modes below a 
certain critical mode while the higher modes are not significantly in- 
volved in the dynamics. Notice, moreover, that an analogous phenomenon, 
i.e. a dynamical involvement restricted to the modes of sufficiently low 
frequency, is know to occur also when a system is coupled to an external 
body, for the case of polyatomic molecules (see [7|). So, it is known that, 
on a short time scale, the high frequencies modes (above fccr) do not get 
dynamically involved, neither by the internal nonlinarities nor by an inter- 
action with external bodies. In both cases a packet of low frequency modes 
is formed which are in mutual equilibrium, and moreover are active in the 
process of thermalization with external bodies. Only after a much larger 
time scale there follows a relaxation of the Fermi-Pasta-Ulam system to 
the true equilibrium state. Before, it appears as if there existed an adi- 
abatic partition (of a dynamical nature) between low and high frequency 
modes. 

If this is the correct picture, it is clear to what property of the Fermi- 
Pasta-Ulam system should our definition of temperature correspond: namely, 
to the mean energy of each of the low frequency modes (those below /ccr)- 
In the rest of the paper we will illustrate the results of some numerical 
computations, which, in our opinion, strongly support the fact that such 
a metaequilibrium is transitive, and that the empirical temperature thus 
defined coincides with the energy of the low frequency modes. 
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Figure 1: Specific harmonic energy of the FPU system, and twice the 
specific energies of the gases, versus time, at high temperature. 

3 The model. 

As mentioned above, the model is constituted by a Fermi— Pasta-Ulam 
system in contact with two bodies, acting, the one as a thermometer and 
the other one as a heat bath. We make the simplest choice, in which the 
two bodies are perfect gases. Each gas is modeled as a system of point 
particles having no interactions among them, while interacting with the 
Fermi-Pasta-Ulam system through some smooth force between each of the 
gas particles and one of the edge Fermi-Pasta-Ulam particles. In more 
detail, concerning the Fermi-Pasta-Ulam system we denote as usual by 
Xi, i = 1,...,N, the distance of the i-th particle from its equilibrium 
position, by Pi its conjugate momentum, and consider the familiar "/?- 
model" Hamiltonian 

1=1 j=0 

involving two positive parameters /3 and fl, with fixed boundary conditions 
xo = 0, xpf+i = 0. Now, to the "left" of the Fermi-Pasta-Ulam chain we 
place a perfect gas which acts as a heat bath: denoting by yi, ttj the i-th 
gas particle's position and momentum respectively, we have —L <yi < xi, 
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with L > playing the role of the volume of the gas. The motion of each 
particle is thus free apart from the fact that it suffers an clastic reflection as 
Ui = —L, and that it moreover interacts with the first Fermi-Pasta-Ulam 
particle xi through a short range potential, which we choose as 

p-iyi-xi)/lo 

(Vi - xi)/lo 

Iq and Vq denoting its range and strength respectively. In agreement with 
the boimd given above, due to the singularity of the potential at yi ~ xi, 
the solutions yi(t) of the equations of motion cannot cross the point xi(t), 
i.e. for all times t one has —L < yi{t) < xi{t). The Hamiltonian of the 
heat bath is thus 



N 2 



,-{yi-xi)/lo 
iVi -xi) 



(2) 



supplemented by the boundary condition that the particles are reflected at 
Vi = L. 
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Figure 2: Same as Figure 1, at low temperature. 

To the ' 'right" of the Fermi-Pasta-Ulam chain we place the thermome- 
ter, which is taken again as a perfect gas, with Hamiltonian 

JV 2 f,-{yi-XN)/lo 

^ 2m [yi -xn) 
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(plus a reflection condition at yi = L) where yi, Wi are the positions and 

momenta of the gas particles respectively; each of the particles interacts 
only with the last particle a; at of the Fermi-Pasta-Ulam system via the 
same potential as for the heath bath. 
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Figure 3: Energy spectrum of the FPU system, at high temperature. 

In our simulations we chose an equal number of particles for the three 
systems, while in principle the number of bath particles should be larger 
than that of the Fermi-Pasta-Ulam system, and this in turn should be 
larger than that of the thermometer. Our choice is dictated only by the 
computational power available: we cannot take the total number of parti- 
cles too large, but at the same time the number of particles in each system 
cannot be too small if a good statistics has to be insured. Taking the same 
number N = 100 of particles for each of the three subsystems, seemed to 
us a good compromise. 

We took as units of mass, length and energy the values m, Iq and 
Vb, which were thus put equal to one in our computations. The values 
of the parameters fl and /3 were set equal to = 400 and /3 = 3742 
respectively. Such strange values come from the following consideration: 
the intermolecular interaction in a crystal is well represented by a Lennard— 
Jones potential, whose relevant parameters (the range and the strength) 
arc of order one with our choice of units. On the other hand the Fermi- 
Pasta-Ulam potential should just be a Taylor expansion of the Lennard- 
Jones potential around the equilibrium position. Performing such a Taylor 
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expansion and putting the parameters equal to one, the indicated values 
for and /3 are found. 

Finally in our numerical simulations we took L — 25; this in order to 
ensure a sufficient total number of collisions (of the order 10^ in our actual 
integrations), while at the same time letting the gas particles be free for a 
large part of their paths. 

4 Numerical results. 

The integration step was taken equal to a twentieth of the shortest period 
Tf = 7r/r2 of the Fermi-Pasta-Ulam chain, and the numerical solutions were 
computed up to times of order 2 • lO^Tf . 

The numerical experiments were performed in the following way. For 
the bath we chose a temperature Ti and took random initial conditions 
extracted from a Maxwellian at the chosen temperature Ti (we also checked 
that the value of the mean kinetic energy should not deviate too much 
from the expected one, in order to avoid too large fluctuations); for the 
Fermi-Pasta-Ulam system we chose initial data at equipartition with a 
temperature Ti/lO and random phases; and finally for the thermometer 
we chose initial data in the same way as for the heath bath, but at a 
temperature Ti/10. We let the system evolve for a time lO^rf, and then 
began to compute the time averages of the kinetic energies of the gases and 
of the harmonic energy of the Fermi-Pasta-Ulam system. The results of 
the computations for two representative cases are shown in Figures ^ and 
El where we report, versus time, the temperatures (i.e. twice the kinetic 
energies per particle) of the gases and the harmonic energy per particle of 
the Fermi-Pasta-Ulam system (actually, time averages of such quantities 
are reported). Let us recall that such three quantities should be equal 
according to the equipartition principle, i.e. for sufficiently long times. 

In Figure n we started with a temperature Ti = 1. One sees that, 
after a time of order lO^rt, the temperatures of two gases and the specific 
harmonic energy of the Fermi-Pasta-Ulam system have become essentially 
equal, although still presenting significant fluctuations. This case should 
correspond to a situation of thermal equilibrium. 

Figure |3 refers instead to the choice of Ti = 0.4. One sees that the 
heat bath and the thermometer still reach the same temperature (although 
after a time almost one order of magnitude larger than before), while the 
specific harmonic energy per particle of the Fermi-Pasta-Ulam system re- 
mains well below the common temperature of the two gases, up to the 
observation time. Actually the curve is so flat that the Fermi-Pasta-Ulam 
system can be expected to possibly reach the equilibrium only on a totally 
different time scale. The global system seems indeed to be in a situation 
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Figure 4: Same as Figure 3, at low temperature. 



of metaequilibrium. 

It appears, however, that the zeroth law can still be valid, and that 
the temperature "measured" by the thermometer is a good empirical one, 
because the temperatures of the two gases have become equal. To under- 
stand to which quantity of the Fermi-Pasta-Ulam system does this mea- 
sured temperature correspond, in Figures 13 and 01 we report the spectra 
(time-averaged energies of the modes versus mode number) of the Fermi- 
Pasta-Ulam system at the end of the two runs. Figure 13 refers to the case 
of complete thermalization, and correspondingly a complete equipartition 
among the modes is found, as expected. More interesting is Figure^ here 
equipartition obtains only among modes of sufficiently low frequency, say 
below kcr = 10, while the energy starts decreasing for larger values of k, go- 
ing down, say for k > 25, to the initial equipartition value 0.04. It does not 
appear as a surprise to observe that the mean energy of the low frequency 
modes essentially agrees with the common temperature of the two gases. 
This seems to indicate that, for our metastable state, the "good" definition 
of temperature of the Fermi-Pasta-Ulam system is the mean energy of the 
(sufhciently) low frequency modes. 
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5 Conclusions. 



In conclusion, we hope to have shown, through our numerical study of a 
Fermi-Pasta-Ulam system in contact with two gases, that there are cases 
of metastable equilibrium for which a notion of temperature can be de- 
fined. However, at variance with the familiar case of equilibrium, such a 
temperature does not coincide with the "canonical" one, namely twice the 
mean kinetic energy per particle. 

As a further comment, we would like to add that the metaequilibrium 
states met in Fermi-Pasta-Ulam systems present characteristics which are 
somehow opposite to those of glasses. Indeed, in the latter case the lack of 
thermalization is ascribed to the low frequency modes and correspondingly 
the thermometer measures the mean energy of the high frequency modes, 
which are the one being in mutual equipartition. 

Finally, we would like to mention that the possibility of having a ther- 
modynamics for situations of metaequilibrium, typically involving the pres- 
ence of adiabatic invariants, was amply discussed in the second part of a 
very interesting paper of Poincare 9^ , which appears to have been almost 
completely forgotten. In fact we became aware of such work only after 
completing the present work, through a conference of V. Kozlov ^U]. In 
fact, V. Kozlov was addressing only the problem dealt with in the first 
part of the paper of Poincare, namely how it occurs that the fast variables 
of an integrable hamiltonian system approach equilibrium, notwithstand- 
ing the reversibility and the return property of the system. In the second 
part of his paper, Poincare was instead considering a situation in which 
one has at first a quick relaxation to a "provisional equilibrium" while a 
"definitive equilibrium" would be attained after a much larger time, i.e. 
one is concerned, in his very words, with "very long times of first order" 
and "very long times of second order" , which is a situation analogous to 
the one discussed in the present paper. 
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